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We calculate the third-order moments of scalar curvature perturbations in configuration space for 
different inflationary models. We develop a robust numerical technique to compute the bispectrum 
for different models that have some features in the inflationary potential. From the bispectrum we 
evaluate moments analytically in the slow-roll regime while we devise a numerical mechanism to 
calculated these moments for non-slow-roll single-field inflationary models with a standard kinetic 
term that are minimally coupled to gravity. With the help of these third-order moments one can 
directly predict many non-Gaussian and geometrical measures of cosmic microwave background 
distributions in the conhguration space. Thus, we devise a framework to calculate different third- 
order moments and geometrical measures, e.g. Minkowski functionals or the skeleton statistic, 
generated by different single-field models of inflation. 


I. INTRODUCTION 

Recent results from ground- and balloon-based exper¬ 
iments as well as from WMAP and the Planck satellite 
have described the temperature anisotropies in the cos¬ 
mic microwave background(CMB) with high precision[T]- 
[S]. This primordial CMB is close to being Gaussian how¬ 
ever there are some features and anomalies that are unex¬ 
plained as shown in Ref. [7] and references therein. The 
study of non-Gaussian contributions in the cosmological 
perturbations is an ideal tool to study the inflationary 
dynamics. Among the most direct data sets for studies 
of inflation and non-Gaussianity are the GMB maps[5HS]. 

Inflation is the initial accelerated expansion of the early 
Universe. Inflationary expansion by more than 65 e- 
folds is needed to explain the observed homogeneity and 
isotropy of the Universe. It is widely accepted that infla¬ 
tion is driven by the potential energy of a scalar inflaton 
field slowly rolling down the potential. Inflation also suc¬ 
cessfully describes the creation of small inhomogeneities, 
needed to seed observed structure in the Universe, as hav¬ 
ing been generated by quantum fluctuations in the infla¬ 
ton field. The quantum fluctuations also got stretched 
and became imprinted on GMB maps and other observ¬ 
ables on cosmological scales. Thus, inflation is able to 
explain not only why the Universe is so homogeneous 
and isotropic but also the origin of the structures in the 
Universe [TOHiB] . Alternatives to inflation have been pro¬ 
posed but no other scenario is as simple and elegant as 
inflation produced by scalar field(s) [T7HT5] . 

In single-field inflationary models generated inhomo¬ 
geneities are described via a single scalar perturbation 
field C (x). The statistical properties of (f (x) are the main 
observable signatures to distinguish inflationary models. 
We know that if the perturbations are exactly Gaussian 
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then all odd n-point correlations functions vanish while 
all even n-point functions are related to the two-point 
function. Thus, in the momentum space the Gaussian 
field C is completely described by the power spectrum 

Pdk), 

(CkCk') = (27r)^<53(k + k')^Pc(fc)- (1) 

Inflation generally predicts that the power spectrum is 
nearly flat P(;{k) oc « const [T51 - [T5] . The fact that 

the observed scalar spectral index n^ = 0.9603 ± 0.0073 
[5] is close to but not exactly unity is also considered by 
many to support the existence of inflation in the early 
Universe. However, many different kinds of inflationary 
models can be made compatible with observations of the 
power spectrum. Thus the study of the non-Gaussian 
signatures is important to reduce the degeneracy in in¬ 
flationary models. Such signatures are contained in non¬ 
trivial higher-order correlations starting with the cubic 
ones. 

Similar to the power spectrum, for three-point correla¬ 
tions one can calculate the bispectrum BQ{ki,k 2 ,k^) as 
a measure of the non-Gaussianity of the initial perturba¬ 
tions 

(CkiCksCka) = (27r)3(5^(ki-bk2-l-k3)H(;(fci,A:2,fc3) . (2) 

The bispectrum carries much more information than 
the power spectrum as it contains three different length 
scales. It is thought that for basic single-field slow-roll 
inflation with a standard kinetic term the non-Gaussian 
effects are small [501, there are many models of inflation 
that give large, potentially detectable, non-Gaussianity. 
Comparing non-Gaussian predictions with observations 
will help us to constrain or rule out different inflationary 
models and give us more insight about the physics of the 
early Universe. 

In this paper we will focus on the study of non- 
Gaussianity through geometrical measures that describe 
visual properties of the initial perturbations viewed as a 
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random field. Examples of standard measures of random 
fields are Minkowski functionals, extrema statistics and 
also more novel measures such as skeleton statistics. The 
simplest Minkowski functional is the Euler characteristic 
or genus density of excursion sets as a function of thresh¬ 
old. It has been shown that for mildly non-Gaussian field 
such local geometrical characteristics can be expressed as 
a series of higher order moments to the perturbation field 
and its derivatives In particular for the Euler 

characteristic [2Tll2^ 


X3d{’^) 
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In the above expression cr^ = (C^)) o"? = ((^C)^), is 
the threshold in units of a and Hi{v) are Hermite poly¬ 
nomials. The first term in the expansion denotes the 
Gaussian part that is proportional to H 2 {v) while other 
terms represent the first non-Gaussian correction in a. 

In this paper we have develop a robust mechanism to 
compute the third-order moments such as (C^)j (C^A^) 
and ((VC)^AC) for single-field models of inflation. This 
links the non-Gaussianity generated by inflation to the 
geometrical observables such as Euler characteristics and 
Minkowski functionals. 

This paper is organized into six sections. In Sec. II, 
we review the theoretical framework of inflationary cos¬ 
mology whereafter we will describe the calculation of the 
three-point correlation function in momentum space and 
the calculation of different third-order moments in con¬ 
figuration space. In Sec. Ill, we will present our numer¬ 
ical technique for the calculation of three-point function 
and briefly discuss different single-field models of infla¬ 
tion with some features in the inflationary potential. In 
Sec. IV, we will present the calculation of moments in 
configuration space while in Sec. V we present the geo¬ 
metrical Minkowski functionals. In the last section we 
will summarize our results and conclude with plans for 
future work. 


II. THEORETICAL FRAMEWORK 


Kline-Gordon equation for inflaton field are given by 

= + Vi4.)) 

4> + 3H,p + V^ = 0 . (5) 

One can define the following slow roll parameters^ and 
the corresponding slow roll conditions as 


e 


H 

H 




(f>) 


These slow roll conditions e <C 1, ^ 1 ensure that the 

inflaton field rolls slowly down the potential and the Uni¬ 
verse inflates for significantly long period. These slow roll 
parameters depend on potential of the inflaton field and 
the model of inflation. For standard single-field inflation 
with quadratic potential V{(j)) = \rri^(fP these slow roll 
parameters are of order 0(0.01) for inflaton field values 
(j) > lOMp. 


A. Calculation of 2-point and 3-point Functions in 
Momentum Space 

In this sections we will present the steps laid down 
by Maldacena to calculate the two-point and the three- 
point correlation function of the scalar perturbations [20] . 
Firstly, one writes the action for the inflaton field given 
in Eq. lousing the Arnowitt-Deser-Misner(ADM) formal¬ 
ism. Secondly, one expands the action to second order 
in perturbation theory for calculation of two-point func¬ 
tion and to third-order for the calculation of three-point 
function. Thirdly, one quantizes the perturbations and 
imposes canonical commutation relations. Next, one can 
define the vacuum state by matching the mode function 
to Minkowski vacuum when the mode is deep inside the 
horizon that fixes the mode function completely. Follow¬ 
ing these steps one can find the power spectrum and the 
Bi-spectrum for scalar perturbations uniiizi- 

In the ADM formalism the space-time is sliced into 
three-dimensional hyper-surfaces E, with three metric 
gij , at constant time. The line element of the space-time 
is given by 

= -N^df + gijidx^ + NUt){dx* + N*dt) (7) 


single-field inflation driven by a scalar field cj) is de¬ 
scribed by the following action in units of = SttG = 

1, c = h = 1) 

S = j d*Xyf^ ( 4 ) 

where V{(j)) is the potential for the inflaton field. In 
Friedmann cosmology with homogeneous and isotropic 
background, the Friedmann equation for scale factor and 


where N and V® and lapse and shift functions. In single- 
field inflation, we only have one physically independent 
scalar perturbation. Thus, we perturb the metric and 
matter part of the action and use the gauge freedom to 


* Our definition of r; = ^ = 2e — 2r;H is commonly used in the 
studies of non-Gaussianity whereas rjn = —ip/{(j>H) is the Hub¬ 
ble slow roll parameter used more commonly in the studies of 
inflation. 
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choose the comoving gauge for the dynamical fields (j) and 
9ij 

S(j) = 0, gij = a2(e2‘^(5y + Uj) , (8) 

where C is the comoving curvature perturbation at con¬ 
stant density hyper-surface 6(j) = 0. In this gauge, the 
inflaton field is unperturbed and all scalar degrees of free¬ 
dom are parameterized by the metric fluctuations (^{t, x) 
while the tensor perturbations are parameterized by tij, 
that is both traceless and orthogonal ditij = t] = 0. The 
conditions in Eq. fixes the gauge completely at non 
zero momentum j20j . The shift and lapse functions are 
not dynamical variables in ADM formalism hence they 
can be derived from constraint equations in terms of 
We shall only study scalar perturbations in this paper. 

Linear perturbation results are obtained if one expand 
the action to second order in perturbation field C 

'S'(2) = J d'^xa^e (9) 

which gives the following equation of motion for scalar 
perturbations v = z( is Fourier space 

Uk -f y ^ Uk = 0. (10) 

where 2 ; = a(p/H and momentum k is in reduced Planck 
units Mpi. This is known as the Mukhanov equation for 
scalar perturbations pH] . Now, one can calculate the 2- 
point function and the power spectrum 

(CkCk') = (2^)3<53(k-k')^Pc(fc) 

Pcik) = ^\nk\^ (11) 

where Uk = Vklz are the Fourier coefficients of C,{x), the 
curvature perturbations. 


To obtain next order results in perturbation theory and 
to calculate non-Gaussianity, one expands the action to 
third-order in scalar perturbations in comoving gauge [201 
[20] • After several integrations by parts and dropping the 
total derivatives one finds the following third-order action 
is often quoted in the literature [20l l29] 

^(3,) = J d'^x + ae^CidCf - 2aeC{d,C){dix) 

+ + ^<d.CMx)d\ + ^e(d^C)id\r 

+ 2/(0^) (12) 

where the last term is variation of quadratic action is 
given by 

■^ = a (“)))“ + X - 5 X = eaX 

/(C) = + terms with derivatives of C- (13) 

Now, the calculation of the three-point function from 
the above action involves integration over the time vari¬ 
able. But, if we use the action in Eq. [^ to calculate 
the three-point function, the integral over time does not 
converge to the end of inflation as pointed out by [30] . 
However, it was shown in[3DH32] that this action can be 
converted into an equivalent form, that gives a conver¬ 
gent three-point function, by adding a total derivative 
term that is given by 

which brings the action to the following form. 


^(3) = j d^x[ a^e(e - + aeX{dCY - - 2aeC(5,C)(9.x) + ^e9^x(a,C)(5,x) + j-^YdX){d.xY (15) 


+ 25(0^^^, <7(0 ^^/^+4 q2^2 


-{dCY + d-Yd^dYdYdjO)] 


1 


Aa?H 


[-{dO{dx) + d-Yd,dYdYd,x))\ 


Here the last term can be eliminated with a field redef¬ 
inition C Cn + 5 (C) because g{C) is only a function 
of derivatives of scalar perturbations C,{t,x) that vanish 
outside the horizon. The above third-order action is an 
exact result without any slow roll approximations thus it 
is even valid for models that deviate from slow roll con¬ 
ditions. Another feature of this action is that it contains 
only first two slow roll parameters e and g while it is 


independent of derivative terms such as 77 '. 

Finally, to calculate the 3-point function in momentum 
space we move to the interaction picture and write the 
Hamiltonian for the action in Eq. [^ as 

P(C)=Po(C)+Pm*(C) (16) 

where Hq is the quadratic part of the Hamiltonian while 
Hint represents all higher order terms in perturbation 
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theory[2D]. The three-point function in the interaction 
vacuum at some time r near the end of inflation is given 
by 


(Cki(T)Ck2('r)Ck3('r)) = 

-i J dr'a([Cki(r)Ck2('r)Ck3(T),-ff*nt(T')]) (17) 


To calculate 3-point function, the interaction Hamilto¬ 
nian Hint is just equal to —5(3) without the integral over 
time coordinate since the conjugate momenta vanish in 
the ADM formalism. Next, we quantize the perturba¬ 
tion held C(x) and dehne the vacuum state. For this 
we expand the Ci^) field into creation and annihilation 
operators and use commutation relations of scalar held 




= 6kk' to get the following result 


(Ck,Ck3Ck3) = *(27r)3 


2 , 2 „, 


dr I — 2 a e M 3 U 2 ^3 




-|-2a^e(e — r])u’[u2u'^ — a^e(2eki.k2 -I- Tyfcgjttj)7/2^3 


,/*,2 k2-k3 


I ^ C) * 1 Z 

U 1 U 2 ki 


3 * /:4c 

-e U1U2 U3 


^2^3 


ki.k2 

kj 


+ c.c. 


-I- distinct permutations n UiiTend) d (^^kj) (18) 


The choice of the vacuum is specihed by the choice of 
mode function Uk selection. This is the main formula for 
the three-point function. This expression will be used in 
the sections to come for the exact numerical calculation 
of three-point function in momentum space. 


B. Non-Gaussian parameters and /ml 


The integral relation for 3-point function given in 
Eq.[^can be analytically evaluated in the slow roll limit. 
Ignoring terms in Eq. 18 gives us the following result 
that was hrst derived by Maldacena|20|. 


(pC\2 

(20) 

i>j 


(Ck3Ck3Ck3) = {27Tys\k, 


A = 


r]- - e 




E 


e 

K 


where * denotes the e and r] values at horizon cross¬ 
ing. The quantity A is a convenient measure of non- 
Gaussianity in the perturbation held. The relationship 
between A and bispectrum is given by 


Byki,k2,k3) 


(27r) 


U.k! 


A. 


( 21 ) 


The bispectrum and A are general measures of non- 
Gaussianity however both these quantities are highly 
scale dependent. Thus, the three-point correlation is 


often described in terms of a local dimensionless non¬ 
linearity parameter /nl- This non-linearity parameter 
/nl was hrst introduced as a measure of local non- 
Gaussianity described by 

C(a^) = Cl{x) - ^/nl (CL(a^)^ - (CL(a^)^)) , ( 22 ) 


where factor I is a matter of convention. For local non- 
Gaussianity, the quantity A is expressed via the non lin¬ 
earity parameter /nl as 


Alocal = ~ ^ ^ ■ 


(23) 


Beyond local model we can dehne a generalized /nl for 
general kind of non-Gaussianity by the following equa¬ 
tion, that also has the advantage of being nearly scale 
independent. 


10Hg(fci,fc2,fc3)n» kf 


(24) 


In this paper we will be using this generalized dehnition 
of /nl instead of A or bispectrum that are both scale 
dependent quantities. The /nl depends on the shapes 
of three-point function triangles and it will depend upon 
scale of triangles if there are features in the inhaton po¬ 
tential. Over the recent years /nl has become a widely 
used measure of non-Gaussianity [33]. 


C. Introduction to Moments 


Non-Gaussianity can also be studied through the 
higher order moments of the perturbation held in con- 
hguration space. Analysis of these moments provide a 
robust measures of non-Gaussianity and has also become 
an important held of investigation [341 135] . In this paper 
we will study how different inhationary models can pre¬ 
dict different non-Gaussian and other geometrical param¬ 
eters in conhguration space. These moments can provide 
important information on the geometrical properties of 
the physical helds, e.g. CMB temperature huctuations, 
and give us non-Gaussianity observables such as extrema 
counts, genus and skeleton]^. 

To calculate the third-order moments we have to take 
the inverse Fourier transform of three-point function in 
momentum space. These third-order moments also con¬ 
tain derivatives of the perturbation held dy and A^. 
The complete and independent set of moments that are 
needed to calculate the observables such Euler character- 
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istic are given by 

<C»M> = / 3 (C-. <:-.<>.) ^■(-■^'■■*'■■1 - (25) 

«^(x)A<(x)) = 1 

xfc2g*(ki+k.+k3).x(26) 

{(VC(i))^A<M> = j 

.(27) 

These moments are always scaled by the correspond¬ 
ing variances a and cti in all physical observables like 
Minkowski functionals and Euler characteristics. 

The details of calculating these moments for different 


inflationary models will be presented in section IV How¬ 
ever, one can the calculate these moments for local non- 
Gaussianity given by Eq. in configuration space and 
there values are given by 




local 


(C^(x)AC(a;)) 


local 


t2^2 


((VC(x))2AC(x)) 


local 


f/NL 

(28) 

f/NL 

(29) 

= 5/nl 

(30) 


The above values of moments are well known results for 
local non-Gaussianity that is independent of any model 
related study of inflation. 


III. NUMERICAL TECHNIQUE 

A. Calculation of Mode Function and Power 
Spectrum 


non-Gaussianity and power spectrum. Thus, we convert 
the Mukhanov equation to perturbation equation for Uk 


u'l = —k^Uk — 2— u' 


(32) 


We solved these equations for background and scalar per¬ 
turbations using the Runge-Kutta method of order four. 
The initial condition for solving these differential equa¬ 
tions are given by the following equations 

a(To) = 1, (j){To) = (j)o (33) 

Mto) = v'kiTo) = ( 34 ) 

We chose initial value of inflaton field (j)Q such that Uni¬ 
verse expands for 70 e-folds for quadratic potential. The 
mode functions originate deep inside the horizon that 
correspond to, our choice of vacuum, the Bunch-Davies 
vacuum. 

The above equations of motion are for single-field infla¬ 
tion with standard kinetic term with any potential V{4>). 
However in this paper we will specifically study quadratic 
inflation V{(j)) = as our base model to check our 

numerics. We will also study two other models that have 
features added to the quadratic model that are known as 
the step and resonance models, details of which are given 
in the next section. 

In Fig we have presented the calculation of dimen¬ 
sionless power spectrum according to Eq. [m The power 
spectrum is mildly dependent on k for the quadratic po¬ 
tential with = — 2e* — T]*. On the other hand for 

the step potential, due to the breaking of the slow roll 
condition because of a sharp step in the potential, we see 
oscillating power spectrum near the step but as we move 
away from the step it follows the same behaviour of the 
quadratic potential (Fig. [^. 

B. three-point function and Cesaro Sum 


To calculate the power spectrum of scalar perturba¬ 
tions we need to solve the background equations of mo¬ 
tion Eqs. [^and the Mukhanov equation Eq. [T0| written 
in terms of conformal time 

a" = 

a 

v'l = — ^ Vk, z = aV^ (31) 

where primes ' donate derivatives with respect to confor¬ 
mal time. These are coupled differential equation with 
first two representing the background and last equation 
for scalar perturbations. Numerically, it is more conve¬ 
nient to work out the differential equation for Uk rather 
then Vk since we finally require Uk = Vkjz to calculate the 


After numerically solving the background equations 
and the equation for scalar perturbations, we insert these 
solutions back into Eq. to calculate the three-point 
function in momentum space. The three-point corre¬ 
lation function is numerically challenging task as it in¬ 
volves integrations that arise from equation (18). The 
integrands consist of three factors of Uk or v!^. multiplied 
by the background factors of a, e and rj. The scalar per¬ 
turbation function Uk oscillates before horizon crossing 
at To, while after horizon it freezes out. Thus, the in¬ 
tegration consists of two parts, before horizon crossing 
(BHC) part and after horizon crossing (AHC) part 


drllr) = J drlir) + 


drlij) (35) 


where tq is the horizon crossing point of the largest k 
mode in the three-point correlation function and /(r) 
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FIG. 1: Power spectrum of scalar perturbations P<;(fc) against 
momenta k for quadratic potential and step potential with 
c = 0.002, d = 0.02 Mpi and 0s = 15.86 Mpi. 


is the integrand of the 3-point function given in Eq. 
that contains background factors and product of three 
oscillating mode functions. The BHC and AHC parts 
of integration present different numerical challenges as 
the first has growing oscillations, as r approaches nega¬ 
tive infinity, while for the AHC part we have to regular¬ 
ize the three-point function by adding a total derivative 
term(Eq. 14) in the action. Without adding this term in 
the action, the AHC part of the integral is divergent as 
one of the term in the initial action grows as the 

scale factor 


The contribution to the integral that arises from before 
horizon crossing poses significant technical challenges. In 
conformal time the initial big bag singularity is pushed 
back in conformal time to r —>■ —oo. Thus, the scalar per¬ 
turbations start deep inside the horizon and keep oscil¬ 
lating till horizon crossing point tq of the largest k mode 
in the 3-point function. Now, there are different meth¬ 
ods to numerically evaluate an oscillating integral over an 
inhnite range. If we cutoff this infinite integral to some 
finite value, due to large oscillations this induces an spuri¬ 
ous contribution of 0(1). Numerically it was shown that 
these kind of integrals can be evaluated by introducing 
an arbitrary damping factor into the integrand but this 
damping factor needs to be chosen carefully [35]. Other 
techniques, such as boundary regularizaion, for evaluat¬ 
ing such integrals are even more complex (SUEZ]. 


We have developed a different numerical technique, 
which is numerically more robust and elegant, using the 
Cesaro resummation of improper series. For oscillating 
integrand /(r), the following expressions gives the defi¬ 
nition of Cesaro integration 


dr/(T) 


lim f dr' 

-^-oo To - r 



dT"l{T'') 


(36) 


This gives a specific definition to the improper integral 
on the left hand side whereas the right hand side is an 
average over the partial integrals that give convergent re¬ 
sult for a wide range of improper integrals [35]. However, 
we extended this method further and we defined a higher 
order Cesaro integral, with one additional average, to 
further improve the convergence as 


lim — f dr' -- 

-^-oo To - T Jpg To - t' 




dT"'l{T"') 


(37) 


In our numerical program we have used this Cesaro in¬ 
tegral with additional average Eq. (371 for faster conver¬ 
gence. 



FIG. 2: three-point correlation function, or generalized /nl, 
integral plotted against the number of e-folds N = ln(a) for 
Equilateral and Squeezed triangles with horizon crossings at 
6.5 and 8 e-folds. The generalized /nl is just the sum of 
the two asymptote(plateau) values at small and large e-folds 
corresponding to the two integrals in Eq|35| 


This method quickly gives convergent results without 
introducing any artificial damping factors. This can be 
seen in Fig. [^ which shows the three-point function in¬ 
tegral result plotted against the number of e-folds for 
an equilateral triangle and squeezed triangle cases. In 
this figure horizon crossings occur tq that corresponds to 
e-folds values of 6.5 and 8 for equilateral and squeezed 
triangle. This Fig [^ describe two different integration 
regimes BHC t < tq and AHC t > tq. In BHC 
regime, we integrate in backward direction from the hori¬ 
zon crossing points using the Cesaro Integral Eq. |37| 
Our technique converges very quickly as can be seen that 
the integral plateaus as we go 5-6 e-folds before horizon 
crossing points. In the AHC regime, we integrate in the 
forward direction that also plateaus soon after horizon 
crossing. Thus, the three-point function integral, or gen¬ 
eralized /nl, is just the sum the two asymptote(plateau) 
values in the before and after horizon crossing regimes 
for each kind of triangle. 
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FIG. 3: Numerical calculation of /nl for an equilateral trian¬ 
gle case(fc — ki = k 2 = fcs) alongside analytical result (Fig. 
a) for quadratic potential. Second plot show results for a 
squeezed triangle fci = 0.005 Mpi while k = k 2 = ks are on 
the x-axis (Fig. b) for quadratic potential. 


To test our procedure, we have calculated the three- 
point function numerically for potential and com¬ 

pared it with the corresponding analytical results given 
by Eq. |20| The numerical results when compared the 
analytical results are plotted in Fig. [^a, b) for equi¬ 
lateral and squeezed triangles. As can be seen that our 
numerics follows very closely the analytical results for the 
quadratic potential. The Fig. |^shows that our numerical 
technique is accurate to below one percent error for equi¬ 
lateral and squeezed triangles. In the extreme squeezed 
triangle limit fc « 5 in Fig. the errors are somewhat 
large as a compromise was made in the calculation due 
to the time and memory constraints. 



FIG. 4: Percentage error in calculation of /nl for an equi¬ 
lateral triangle case when compared with analytical results 
plotted against fc = fci = A :2 = fcs. While for squeezed tri¬ 
angle percentage error is plotted against k = k 2 = k^ while 
fixed fci = 0.005 Mpi for quadratic potential. 


C. Application to Models with Significant 
Non- Gaussianity 

Single field inflation with quadratic potential, V ((/>) = 
^rn?(jp, is categorized as chaotic inflation. The dynam¬ 
ics of inflation can be solved exactly in the slow roll limit 
for quadratic potential, thus we have used it as our base 
model for our study. However, it is was shown by Mal- 
dacena that single held slow roll models of inhation give 
rise to small non-Gaussianity of 0{e) as can seen in Fig.[^ 
and Eq. In this paper we take the inflaton held mass 
to be TO = 

In single-held models of inhation signihcant non- 
Gaussianity can arise from different kinds of non trivial 
potentials, for example adding some features to quadratic 
potential. One such model consist of localized breaking 
of slow roll conditions. This can be achieved be adding 
a step in the quadratic potential as proposed by [35] 

Vstep{4>) = + ctanh • (38) 

where c is the height of the step and d the width of the 
step centred at (j)s- This model has been used to im¬ 
prove the ht between LCDM and observed power spec¬ 
trum and we have used the best ht parameters for the 
step potential|35|. When inhation rolls down through 
this step it goes through a sharp acceleration in inha¬ 
tionary dynamics. The three-point correlation function 
is proportional to e and 77 parameters and deviation from 
slow roll gives rise to large interaction of modes near the 
step. This gives rise to large non-Gaussianity and gen- 
eralized /nl ~ ^ [36] . It important to mention here 

that e does not change much by the step potential how¬ 
ever f] « gives dominant contribution to the 3-point 

















function. Now, the three-point function integral Eq. 18 


gets most of its contribution near the horizon crossing 
part of the modes. Thus, due to the step in the poten¬ 
tial, the modes that cross the horizon near the step get a 
sharp kick that gives rise to large non-Gaussianitv [571110] . 

In Fig. we present the calculation of /nl for equilat- 
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FIG. 5: /nl for the step potential model plotted against k = 
k\ = k 2 = fcs for equilateral triangle case with step height 
and width given by c = 0.002 and d = 0.02Mpi respectively 
while step is located at (f>s = 15.86 Mpi. 
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FIG. 6: /nl for resonance model plotted against fc = fci = 
k 2 = fcs for equilateral triangle case for c = 5 x 10“^ and 
0.0007. 


eral triangle for the step potential model with parameters 
c = 0.002, d = 0.02 Mpi and t/fg = 15.86 Mpi. The mode 
function Uk that exits horizon at ^ = (/g corresponds 
to ks « 0.01 Mpi. It is noted that peaks in /nl plot 
for the step potential are of 0(10) in amplitude that is 
about 500 times the /nl of quadratic potential. How¬ 
ever, if we move away from the step fcg, fc > 0.5 Mpi and 


k < 0.005 Mpi, the value of /nl for step potential tend 
to approach values for the quadratic potential. 

The other kind of models that can give rise to large 
non-Gaussianity are models with global features in the 
potential like small oscillation on top of quadratic poten¬ 
tial also know as Resonance mo del [37] 

Vres{4>) = (1 + csin((//A)). (39) 

where c is the amplitude of oscillation and A is the fre¬ 
quency. In this model the non-Gaussianity is generated 
in the sub-horizon scales when modes are deep inside 
the horizon and they resonate, or interfere constructively, 
with the physical frequency oj « This model intro¬ 
duces a new scale A into the quadratic potential as there 
small ripples in the entire plane of the potential. These 
ripples in the potential causes oscillation in k space in 
the generalized /nl with amplitude (see Fig. 

while the power spectrum remains almost flat with tiny 
ripples [57]. This kind of mechanism may be realized in 
brane inflation with duality cascade [511142] . Fig. shows 
our computational technique works well for this kind of 
models as well. 


IV. CONFIGURATION SPACE MOMENTS 


Our goal is to develop a numerical technique for calcu¬ 
lating the configuration space moments for perturbation 
/ field in a general single-field model of inflation. For that 
we need to integrate three-point function in momentum 
space as given by Eqs. (^27). As a starting point we 
shall consider slow roll models where these moments can 
be calculated analytically for the case of flat power spec¬ 
trum « const and then use the acquired insight to 
develop a general numerical procedure. 

For instance, to calculate the (/^(x)), we substitute 
Eq. [T^ into Eq. to get 


/ n -f k2 -f- ka 


2=1 


47r^ 


“Tki, 


k2,kf 


(40) 


where A is given by Eq. |20| For near flat spectra such in¬ 
tegrals generally have both infrared and ultraviolet diver¬ 
gences, thus the correspondent moment is formally infi¬ 
nite. Ultraviolet divergences are regularized by the finite 
resolution of our measurements and should be studied in 
the context of a specific experiment. Infrared divergen¬ 
cies, on the other hand, seem of more conceptual nature 
since they come from contribution of the modes much 
larger than the observed Universe, which potentially re¬ 
flect also less understood physics. To study them we 
introduce infrared kmin and ultraviolet kmax cutoffs and 
consider asymptotic behaviour as z = kmin/kmax —>■ 0. 
This procedure is somewhat ambiguous as to the order of 
imposing the cutoff and integrating over the (5-function. 
This ambiguity does not, however, affect the coefficient 
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of the leading infrared divergent term, which is, as we 
will see, what matters. It does affect subleading terms 
and this freedom can be used to minimize them in the 
numerical calculations. We suggest the procedure to use 
the ^-function to eliminate in every term in the integral 
the least infrared divergent momentum and then integrate 
over the angles and the magnitudes of remaining two mo¬ 


menta in kmin and kmax limits. 

Using the above mentioned procedure to evaluate the 
momentum space integrals, we also calculated the other 
two moments ((^(x)A((x)), ^(VC(x))^A^(a;)) by sub¬ 
stituting Eq. [T^ in to Eqs. [^ and The asymptotic 
results are given below up to order 


(C"(^)) 

{C{x)AC{x)} 

((VC(a:))2AC(x)) 


(Pk)" + 2e*) \n{z) - 3eHz) + e (15 - 241n(2)) + 0{z^)^ , 

-{P^kfkl^x [{ri* + 2e*) ln(z) -j{7t^- 28 + 32 ln(2)) + 0{z^)^ , 


(41) 

(42) 

(43) 


Evidently, the first two moments are divergent in z —> 0 
limit, (C^(a;)) oc In^(z) while (^C'^{x)A({x)) oc ln(z). The 
modified regularization sequence retains the coefficients 
in the divergent terms but changes the constant ones. 
Here we note that the last moment with the most deriva¬ 
tives ((VC(a;))^AC(cc)) is finite as z —>■ 0 and thus remains 
sensitive to regularization details. 

Now, if we look at the expression for Euler character¬ 
istic (Eq. or any other observables, the configuration 
space moments are alw ays divided pS] by the variances 
and af given in Eqs. 

a^ = {e{x)} = P^,ln(^) (44) 

\ ^min / 

al = ((VC(a;))2) = . (45) 


Thus, even though the moments by themselves are diver¬ 
gent, the observed quantities are finite. Normalized or 
scaled by these variances the moments to leading order 
in slow roll are given by 


^3 

^3 

U 3 


+ 2 e*), 

{(^{x)AC,{x)) 


t2^2 


~ + 2 e*), 


{{VC,{x)fAC,{x)) 



81n(2) 



(46) 

(47) 

(48) 


For numerical calculation of momentum space integrals 
finite momentum space cutoff is unavoidable. The above 
analytical results guide us to the following numerical pro¬ 
cedure for the general single-field model. We calculate 
the moments using the exact result of 3-point function 
Eq. with the triangularity condition ki -|- k 2 -f k 3 = 0 
applied to the least divergent momentum.To numerically 
obtain z —>■ 0 result we integrate the remaining two mo¬ 
menta in finite range [kmimkmax], then vary and 

find the asymptotic limit (plateau value) for the scaled 
moments T 3 and U^. Procedure we follow is that we 
fix kmin = 0.005 Mpi that corresponds to the mode that 
inflates 64 e-folds after horizon crossing that is roughly 
the largest scales of the observed Universe and vary kmax 
or z. As seen in Fig. for quadratic potential the nu¬ 
merical calculation of S' 3 , and U 3 , gives very stable 
result already by kmax ~ 0.1 Mpi. This fast convergence 
reflects smallness of the constant terms in Eqs. and 
—1.6e and 0.8e correspondingly. Slight drop in values 
of moments in Fig. is a numerical artifact as higher 
values of kmax > 2Mp/ require very fine resolution Afc in 
momenta integrations. 

Thus we see that the infrared divergences for near flat 
spectra can be dealt with very efficiently during numeri¬ 
cal momenta integration. It is sufficient to integrate over 
just a decade of wave-numbers to obtain a good approx¬ 
imation to the asymptotic values. The quadratic poten¬ 
tial gives the values S 3 = 0.0517, T 3 = —0.0664 and 
U 3 = —0.0264 with at most 5% error in the calculation 
of these moments. 


Note that the numerical coefficient in front of e* in U 3 
is also equal to two with better than 2 % accuracy, thus 
all three normalized moments are determined by rj* + 
2e*. From the above moments, we have found that the 
local non-Gaussianity parameter for single-field slow roll 
inflation is /nl = + 2 e*) which is independent of 

the shape and scale of the triangles. 


For potential with features such as step potential, the 
answer for the moments depends on the range of k- 
integration in relation to the modes that give additional 
contribution to non-Gaussian signal. For the step po¬ 
tential, as Fig. demonstrates, the support for correc¬ 
tion to the /nl is finite, encompassing the range from 
k « 0.01 Mpi to fc « 0.2 Mpi for the parameters used 
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FIG. 7: Moments plotted against kmax with fixed kmin = 
0.005 Mpi for quadratic potential. 



^max 


FIG. 8: Moments plotted against kmax with fixed kmin. = 
0.005 Mpi for step potential with c = 0.002, d — 0.02 Mpi and 
0 . = 15.86 Mpi. 


creasing as the correspondent powers of the variances, 
i.e. as Ija^ « ln“^for S' 3 , « 

for T 3 and Ijol « for U-^. 

This explains the quick convergens of T 3 and C /3 mo¬ 
ments to their baseline values as kmax > 1 Mpi and slow 
logarithimic decrease of S 3 exhibited in Fig. 

In cosmological applications the kmin to kmax range 
depends on observational setup or analysis choice. For 
large-scale structure studies the ratio of largest observ¬ 
able scale ^ 10 Gpc to the smallest presently mildly non¬ 
linear scale ^ 10 Mpc is of order of a thousand, which 
points to kmax 5 Mpi. This is the value that we adopt 
in the following sample calculations. 

For the step potential we plot the moments against 
the parameters c and d as displayed in Figs. and 
The limit of the quadratic potential is reached as c —>■ 0 
or d —)■ 00 . It can be seen that the magnitude of S 3 
moment increases significantly as the height of the step 
c increases (Fig[^ or the width d decreases (Fig. 10 1 . At 
the same time, the moments T 3 and U 3 remain practically 
unchanged from their small values in quadratic potential 
limit, as expected from our analysis of kmax behaviour. 
The first moment S 3 is plotted against c and d in Fig. 
which shows that sharper the step, smaller width and 
higher step, the higher the moment. 



0.001 0.01 
c 


there. Fig. shows the change in the moments if we 
fix kmin = 0.005 Mpi but vary kmax- The asymptotic 
behaviour can be understood from the following consid¬ 
eration. Let us write 

. Jc 1 — oQ^a-d ^ (O /qq'i 

*^3 ya-mimi^max) — *^3 1“ rr r, 

^ y^mini^max) 

and correspondingly for the other moments, where 
A signifies additional contribution to the momen¬ 
tum integral above the baseline quadratic values. In 
such a split the first term, is practically inde¬ 

pendent on the kmin or kmax, as we discussed above. 
The second correction term depends on the integration 
range. As soon as this range encompasses the support 
for the extra non-Gaussian contribution, it starts de- 


FIG. 9: Absolute values of the moments for the step potential 
plotted against parameter c for fixed d = 0.01 Mpi. 


V. GEOMETRICAL STATISTICS AND 
MINKOWSKI EUNCTIONALS 

Many geometrical statistics, including Minkowski func¬ 
tionals, of the mildly non-Gaussian fields can be ex¬ 
pressed as series expansion in higher-order moments 
of the field and its derivatives with powers of a con¬ 
trolling the expansion order [23]. First order non- 
Gaussian corrections are linear in a and defined by cu¬ 
bic normalized moments that we have studied earlier. 
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FIG. 10: Absolute values of the moments for the step poten¬ 
tial plotted against parameter d for fixed c = 0.01. 





(C3\ 

FIG. 11: Moment S3 = for the step potential plotted 
against parameters c and d showing that this moment peeks 
when the step is sharpest. 


For the models with nearly flat power spectrum that 
we consider, we adopt the value of the variance a « 
3/Pcln{kjnax/kmin) « 1-2 X lO”'^ with = 2.2 X 10“® 
suggested by Planck data [ 8 ] and kmax/kmin = 10^. 

The simplest geometrical statistic or Minkowski func¬ 
tional is the filling factor fv = dC,P{(), i.e. the frac¬ 

tion of volume above the threshold v. Its moment expan¬ 
sion gives 


fviv) = ^Erfc -h ^^aH 2 {v)Y + Oia^). (50) 

The first non-Gaussian correction only depends on S 3 
multiplied by Hermite polynomial of order two. Thus, 
linear in a non-Gaussian part of fv{^) has a global shape 
that is independent of any model while its magnitude 


will depend on the magnitude of moment S 3 . Thus for 
the step potential, the Ailing factor can be as large as 
500 times the fv for quadratic potential. In Fig. a 
single curve shows the non-Gaussian part of the filling 
factor both for the model with quadratic potential with 
fv correction of order 0 ( 10 “^) (as labeled on the left 
vertical axis) and for the step potential for which it is 
0 ( 10 “®) (the right vertical axis). 



FIG. 12: Non gaussian part of filling factor fv as the func¬ 
tion of threshold u. Values on the left vertical axis are for 
quadratic potential while values on the right axis are for step 
potential with c = 0.01 and d — 0.01 Mpi. 


There are several advantages to use the value of the Ail¬ 
ing factor fv instead of as a variable in which to express 
all other statistics. Indeed, the fraction of volume occu¬ 
pied by a data set is often available even from the limited 
data, whereas specifying v requires prior knowledge of the 
variance cr which may not be easily obtainable. In some 
cases non-Gaussian analysis itself gives more robust way 
to determine the variance. Following |4r{j we introduce 
the effective threshold Vf = •\/2Erfc“^(2/v) be used as an 
observable alternative to v. To the first order correction 
in a we have the relation 

(51) 


Another Minkowski functional in 3D is the area (per 
unit volume) of isodensity contours N 3 {v). Up to first 
non-Gaussian contribution, this quantity is expressed in 
terms of third-order moments as 


M3{v) 


^2 

2 e~^ a I 


1 + cr 



H3{y) + 



,(52) 


depending on moments S 3 and T 3 . In Fig. we show 
the non-Gaussian A /3 (^) for quadratic and step potential. 
Besides vastly different amplitudes, it exhibits different 
shape that distinguishes the two models. For quadratic 
potential H 3 contribution is small and A /3 has just two 
measurable extrema with small secondary ones near the 
edges, while due to the large value of S 3 the H 3 part 
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V 

FIG. 13: Non gaussian part of area of isodensity contours 
Afsiv) for quadratic potential in black with values on left 
vertical axis as a function of threshold v. Similarly, N 3 {v) 
for step potential, for c = 0.01 and d = 0.01 Mpi, is show 
in red(lighter colour) with values on right vertical axis as a 
function of threshold v. 



y, 

FIG. 14: Non gaussian part of area of isodensity contours 
J^ 3 {vf) for quadratic potential in black with values on left 
vertical axis as a function of filling factor Uf. Similarly, 
Af 3 {v>f) for step potential, for c = 0.01 and d = 0.01 Mpi, 
in red(lighter colour) with values on right vertical axis as a 
function of filling factor Uf. 


is prominent for the step potential, giving rise to four 
distinct extrema, as seen in Fig. [T^ 

As a function of i/f, i.e of the filling factor, the area of 
isodensity contours is 


Msivf) 


^2 

2e“T cTi 

'/Sna 


1 — a 





(53) 


which demonstrates the general outcome of eliminating 
the highest order Hermite polynomial term when switch¬ 
ing from 1 / to Vf. Because of that, in Fig. 14 the area 
of isodensity contours has a model independent 


shape that is proportional to e~^Hi{vf). The 
for the step potential has values of order 10“^ on right 
axis while for quadratic potential it has values of order 
10 “® on the left axis. 

An important Minkowski functional in cosmology is 
the Euler characteristic (or genus) that is used to char¬ 
acterize the topology of the isocontours of random fields. 
The definition generally used by cosmologists is that the 
genus is the number of holes minus the number of isolated 
regions above a threshold v in a random field (thus, for 
one isolated region, the genus is just the number of holes 
minus one), while Euler characteristic of the excursion 
set is just minus genus, see Eq. Euler characteristic 
density can be also considered as a function of filling- 
factor-deduced threshold Vf rather than the threshold v. 




1 

w? 



a 








(54) 


We have calculated Euler characteristic xsp of the 
3D perturbation field C for two models of inflation, 


the quadratic inflation and step potential model. For 
quadratic potential, the non-Gaussian part of xsd as a 
function of threshold i/ is plotted in Fig. showing the 



-4 -2 0 2 4 


y 

FIG. 15: Non gaussian part of Euler characteristic X 3 d as 
a function of threshold v for the quadratic potential. The 
shaded region shows the error bars coming from the error in 
calculation of moments that are roughly 5 percent. 


small amplitude of the non-Gaussianity of order 10 ^ 
and the presence of all three Hi, harmonics. Eu¬ 


ler characteristic as a function oi Vf is given in Fig. 16 
For the quadratic inflation 53 -I- 3 T 3/4 « 0.002 is notamy 
smaller than S 3 + 9173/4 « —0.008, hence the result is 
dominated by the Hi term that has only one zero cross¬ 
ing at origin as can be seen in Fig. [THl 

While the 0(10“^) non-Gaussian correction to Eu¬ 
ler characteristic for the quadratic potentialis small, as 
expected, and is hardly observable, for the step poten¬ 
tial, it will have the magnitude 10 to 500 times larger. 
Fig. 17 shows X 3 d(j^) curves for two different sets of 




























13 



-4 -2 0 2 4 

Vf 



FIG. 16: Non gaussian part of genus as a function of filling 
factor Vf for the quadratic potential. The shaded region shows 
the error bars coming from the error in calculation of moments 
that are roughly 5 percent. 


FIG. 18: Non gaussian part of genus as a function of filling 
factor Vf for the step potential for two parameter sets (c = 
0.0018, d = 0.022 Mpi in black on left axis) and (c = 0.01, 
d = 0.01 Mpi in red, lighter colour, on right axis). 


c and d parameters, with magnitude of the effect dif¬ 
fering by an order of magnitude, namely 0(10“®) for 
c = 0.0018 and d = 0.022 Mpi and 0(10“^) for c = 0.01 
and d = 0.01 Mpi. At the same time the shapes of 
the Euler characteristic curves are very similar domi¬ 
nated by SzHfyiyv) term. This can be seen explicitly 



V 


0.000019 


0 


-0.000019 


FIG. 17: Non gaussian part of genus as a function of threshold 
V for the step potential for two parameter sets (c = 0.0018, 
d = 0.022 Mpi in black on left axis) and (c = 0.01, d = 
0.01 Mpi in red, lighter colour, on right axis). 


also when measuring Euler characteristic as a function 
of Vf as shown in Eig. 18 With S 3 dominant for all 
the values of step parameters c and d studied, the non- 
Gaussian part of Euler characteristic exhibits a universal 
form y3D oc ( 7 S 3 (Hi{vf) -\- H 3 (vf)) that is a distinguish¬ 
ing signature of this particular model of non-Gaussianity. 


VI. RESULTS AND DISCUSSION 


The main goal in this paper was to develop theoreti¬ 
cal formalism that links early Universe inflationary mod¬ 
els to the observable geometrical characteristics of the 
initial field of scalar adiabatic cosmological perturba¬ 
tions. The link to non-Gaussian features in such statis¬ 
tics as Minkowski functionals, extrema counts, and skele¬ 
ton properties is provided by studying the higher order 
moments of the perturbation field and its derivatives in 
configuration space. 

We have investigated from the first principles the third- 
order configuration space moments that give first non- 
Gaussian corrections. To calculate these moments we 
have to calculate the three-point bispectrum in momen¬ 
tum space and then integrate over the three momenta 
with appropriate combinations of derivative operators. 

We presented a complete prescription for the calcu¬ 
lation of three-point function, in momentum space, for 
general single-field models of inflation, improving on the 
computational procedures used in the previous literature. 
We have devised a method to precisely calculate the re¬ 
quired time integral over products of perturbation mode 
functions in “before horizon crossing” as well as in “after 
horizon crossing” regimes. Eor that we stressed the neces¬ 
sity to use the action that gives stable three-point func¬ 
tion integral result till the very end of inflation [30] and 
we developed novel technique using modified Cesaro sum¬ 
mation to integrate the highly oscillatory “before horizon 
crossing” part of the three-point function integral. Thus, 
we calculated the bispectrum in terms of generalized /nl 
that even works for models that break slow roll condi¬ 
tions and where there is significant variation in slow roll 
parameters. 

In the configuration space we have studied infrared 
and ultraviolet divergences of the third-order perturba¬ 
tion moments. Analytical analysis of the slow-roll models 
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allowed us to illuminate the role of the dominant infrared 
divergent logarithms in the moment calculations. These 
terms determine observable moments scaled by the ap¬ 
propriate powers of the variance of the field or its deriva¬ 
tive. The moments that do not contain infrared diver¬ 
gence are found to be sensitive to ultraviolet smoothing 
prescription. We implemented this understanding in nu¬ 
merical procedure that transforms bispectrum results to 
the configuration space moments. 

To calculate the moments we had to integrate out tri¬ 
angles of all shapes and sizes. We introduced the reg¬ 
ularization in the momentum space integrals with both 
infrared kmin and ultraviolet k^ax cuts. We observed 
that for the basic flat spectrum these moments do not 
depend on instead leading logarithmic terms de¬ 
pend on the ratio z = of the smallest the largest 

scales and can be robustly select relative to subdominant 
term by varying this ratio which is numerically easier to 
do by varying kmax- It was further noted that these mo¬ 
ments divided by their corresponding variances are finite 
quantities and they reach a plateau for large values of 
kmax hence they are also independent of physics of the 
smallest scales (see Figs. [^and[^. For the step poten¬ 
tial it was shown that the moment S 3 is the one that is 
affected the most while the other two moments T 3 and 
U 3 are very mildly affected by the step if the range of 
k integration encompasses all the modes with enhanced 
non-Gaussian response. 


The Minkowski functionals of the cosmological pertur¬ 
bation field generated on inflation depend on the bal¬ 
ance between the moments of the field and its deriva¬ 
tives up to the second order. We have demonstrated typ¬ 
ical behaviour of the non-Gaussian corrections to several 
Minkowski functionals of ( field in the slow-roll quadratic 
model and the model with step potential. Such correc¬ 
tions show different signatures for different models of in¬ 
flation. In particular step potential, as representative of 
the class of models with enhanced non-Gaussian response 
over limited wavelength range, leads to a characteristic 
universal form of the genus curve as function of the filling 
factor if measurements include a sufficiently large interval 
of scales. 
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